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Abstract 

We derive the energy asymmetry given the proto-neutron star during the time 
when the neutrino sphere is near the surface of the proto-neutron star, using the 
modified URCA process. The electrons produced with the anti-neutrinos are in Landau 
levels due to the strong magnetic field, and this leads to asymmetry in the neutrino 
momentum, and a pulsar kick. Our main prediction is that the large pulsar kicks start 
at about 10 s and last for about 10 s, with the corresponding neutrinos correlated in 
the direction of the magnetic field. 

1 Outline 

Supernova gravitational collapse of massive star often creates a pulsar. 

Many pulsars are moving with velocities much larger than stars in our galaxy: 
PULSAR KICK 

I'll review hydrodynamics of the collapse and neutrino emission < t < 10 s: 
NO GO 

I'll review sterile neutrinos possibly seen at LSND (&: MiniBooNe ?) giving 
the pulsar kicks. 

Research by Ernest M. Henley, Mikkel B. Johnson and Leonard S. Kisslinger: 
We get pulsar kicks consistent with observations from asymmmetry due to the 
electrons produced in the modified URCA process in strong magnetic fields being 
in Landau levels 

2 Review Supernovae and Neutron Star Formation 

A supernova is the gravitational collapse of a massive star (> 8 sun masses): 

1. Collapse to density > 10 14 g cm -3 > nuclear density 

Shocks, bounce, etc. 

2. Protoneutron star formed ~ 0.01 s. 

Neutrinos trapped in neutrinosphere. Radius of neutrinosphere ~ 40 km. 
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Later stage t ~ 10 s 



3. From 0.1 to 10 sec neutrinosphere contracts from ~ 40 km to protostar radius ~ 10 

km 

Neutrinos carry gravitational energy from the emerging star. 

4. From ~ 10 to ~ 50 sec n-n collisions dominate neutrino production and protoneutron 
star cooling. 

The modified URCA process dominates energy emmision by neutrinos. 

n + n — > n + p + e~ + v e 

2.1 Pulsar Kicks 

Resulting neutron stars (pulsar) have a bimodal velocity distribution 

< v < 300km s _1 
1000 < v < 1500km s _1 



(see, e.g., B.M.S. Hansen and E.S. Phinney, |astro-ph/9708071) 



I.e., pulsars with high luminoscity tend to have velocities more than an order of magnitude 
greater than star velocities in our galaxy. 
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2.2 Attempts at Explanation 



There have been many papers written on attempts to explain the pulsar kicks by processes 
taking place during the first 10 s, when most of the energy from the collapse is carried 
off by neutrino emission. The main attempts have been asymmetries in the hydrodynamic 
processes associated with the collapse, the scattering of neutrinos produced by the URCA 
process by nucleons polarized by the strong magnetic field within the neutrinosphere, and 
by sterile neutrino processes. We briefly review these attempts. 

2.2.1 Asymmetries in core collapse 

There have been many attempts to explain the pulsar kicks by the asymmetries in the initial 
collapse of the massive star. Recently authors of a calculation postulating a slow collapse (t 
> Is) claim that they can predict pulsar velocities greater than 500 km/s (L. Scheck et. al., 
RR.L. 92 (2004) 011103). In a review by J. Murphy, A. Burrows, and A Heger, Astrophysics 
J. 615 (2004) 615, a careful analysis of all hydrodynamic calculations shows that one cannot 
obtain a kick of more than about 200 km/s. 

2.2.2 Processes involving standard neutrinos with a strong magnetic field 

During the first 10 s, when the neutrinosphere has a radius starting about 40 km, the URCA 
process dominates neutrino emission (see O.F. Dorofeev et al, Sov. Astron. Lett. 11 (1985)), 

p + e~v 
n + u e 

P+ Ve , 

where the nucleons are polarized by the strong magnetic fields of the protoneutron star. A 
few years ago calculations of elastic neutrino scattering from the polarized nucleons ( C.J. 
Horowitz and J. Piekarewitz, N.P. A640 (1998) 281; C.J. Horowitz and G. Li, PRL 80 (1998) 
3694) found this process to be promising to explain the pulsar kicks. 

However, if one includes the neutrino absorbtion opacities modified by the strong mag- 
netic fields, the asymmetric neutrino emission is reduced, and these processes with standard 
neutrinos cannot account for the pulsar kicks, D. Lai and Y-Z. Qian, |astro- ph / 9802345 , 
ApJ.(1998). Note that although the neutrinosphere at 1 s is about 40 km, the mean free 
path of neutrinos is only or the order of 1 cm, and the asymmetric neutrinos are not emitted. 

Mechanisms with standard neutrinos during the first 10 s have not been able 
to account for the large pulsar velocities. 



n 

e~ +p 
e + + n 
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2.2.3 Pulsar kicks from sterile neutrinos in a strong magnetic field 

If sterile neutrinos exist they can be produced by oscillations in the protoneutron star matter. 
In the presence of the strong magnetic fields of the protoneutron star neutrino asymmetries 
will be produced. Since the opacities of sterile neutrinos are very small, this has been studied 
as a possible source of the pulsar kicks (A Kusenko and G. Segre, PRL 77 (1996) 4872; PL 
B396 (1997) 197. 

Using an extension of usual concept of the CKM mixing matrix giving the weak interac- 
tion states in terms of the mass eigenstates, sterile/active neutrino mixing given by mixing 
angle \9 m \ « 1: 

\u\> = cos6 m \u e > —sind m \v s > trapped 
1 1*2 > = sin6 rn \u e > +cos6 rn \u s > not trapped 

In a calculation that requires the consistency of sterile neutrinos and dark matter (G.A 
Fuller, A. Kusenko, I. Mocioiu and S. Pascoli, PR D68 (2003) 103002) the parameters needed 
to give the kick were found, as illustrated in the figure below 



10 

m s (KeV) 



1 

—7 8 9 

10 10 10 

sin 2 9 

One sees that a sterile neutrino mass of 1- 10 Kev is needed, which is not expected. 

New data on the sterile neutrino's mass and mixing angle is expected from 
MiniBoone. If the results are similar to those of LSND, with fits such as Am 2 = 
0.2ew 2 and sin 2 9 = 0.068 (From William Lewis, LANL) there is a possibility to get 
the pulsar kicks from sterile neutrinos. 
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Fuller et. al. PR D68 (2003) 103002 
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3 Pulsar Kicks From the Modified URCA Process in 
a Strong Magnetic Field 

In the later stage (t >~ 10s) the Modified URCA process: 

n + n — > n + p + e~ + P e (1) 
dominates the cooling of the neutron star. 

URCA DIAGRAMS 

P l P 3 P l P 3 P 2 P 4 



^> — n n > 1 n n > . 5> — p 

(7t°) (7t ); (7C _ ) ; 

k v k y k Y 

P 2 ^\ P 4 P 2 q \ P 4 P, ' P 3 

(1) V e \ ~ C (2) (3) 



( 




The early work includes 

J.N. Bahcall and R.A. Wolf, PRL 14 (1965) 343; PR 140 (1965) B1452; 
E.G. Flowers, P.G. Sutherland and J.R. Bond, PRD 12 (1975) 315. 
B.L. Friman and O.V. Maxwell, Astrophysical J. 232 (1979) 541, who used 
pion exchange and short-range potential in detailed calculations. 
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Modified URCA ingredients: Nonrelativistic OPE interaction and Standard 
W-S Model: 



Vope = (^faU.k—^aW.kTM.TW, (2) 
rrin k z + m~ 

with / ~ 1.0, <j,t the Pauli spin, isospin operators, and (1,2) refer to the two nucleons at 
the pion vertices. 

Weak axial n-p vertex in nonrelativistic Standard Weinberg-Salam model: 

G - 

Weak Axial W A = --^g A Xp • <?Xn , (3) 
with G=^t,9a = 1-26, and the lepton current defined by 

We need the matrix element product 

V I A/f 1 2 = V IT/ Tyf/ 1 2 

'—'lepton, nucleon spins \ lvl \ — / -'lepton,nucleon spins \ v strong' v weak] 

In previous research the vector-axial vector product was used. However, we realize that 
in a strong magnetic field the axial- axial product, which is much larger, can be used, as the 
magnetic field provides the asymmetry. 

We also make use of the fact that in the strong B field with the electrons in Landau 
levels, u(q e ) — > ^ LanAm (g e ) one can get a large pulsar kick without n-polarization, greatly 
simplifying the calculation. The stationary wave function of an electron in a strong magnetic 
field is (J.J. Matese and R.F. O'Connel, Phys. Rev. 180, 1289 (1969)) 

B = Bz 

4>? andau ' n= °(qi,p z ,cf>) = i( y /T)- 1 e-W*' 

iJ Dirac >-(q e )S(q e -p e z ), (4) 

where 7 = Bm 2 e /2B C) with the critical magnetic field B c ~ 4 x 10 13 C I.e., in a strong 
magnetic field the electron moves as a plane wave in the z direction, with energy of Landau 
levels, and is fully negatively polarized in the lowest Landau level. 
The leptonic traces needed are 

TrileptonWllj] = Tr^a) + Trp^-]^) 
+Tr[4/,]( 7 ) 

Tr[/U 7 ,y(a) = -47-^ -^(l + a-fe)) 
TrfeP) = *47-%- fc « 

Trfe](7) = 47-W + ^n- (5) 
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Our notation for the matrix element products is 

k = p x - p 3 , p = p 2 - pi , ;k e = k + p 

A = ( ) —j=9A- 



a/2 uj(k 2 + ml) 
R(k) = f,t m2 \ ■ ( 6 ) 



3.1 Neutrino emissivity 

Defining the Axial Product Matrix Element, Maa = \Ma\ 2 '- 

M aa = iTrilepton^ljljlWAiVoPE 

+exchange)]ij\ 2 , (7) 

the general form of neutrino emissivity with A-A proces is 

II^ 2^(2^)3(2^)3^^ 

— Ei n itial)5(p fi na l ~ P initial) M A A? , (8) 

where T is the product of the initial and final Fermi-Dirac functions corresponding to the 
temperature and density of the medium. J 4^ omitted in present work. 

The nucleons and the electrons are in thermal equilibrium. A crucial observation of this 
work is that the time scale of the strong interaction in the modified URCA process is short 
compared to the time scale for the electron, which is that of electromagnetic interaction. 
Therefore the proton quickly reaches thermal equilibrium, while the process of the electron 
state transforming to the n=o Landau state does not intefere with the proton reaching its 
Fermi momentum. Therefore we can use the values for the nucleon and electron momenta 
in the matter of the protoneutron star, as derived in the Friman-Maxwell paper. 

The J d 3 q e , involving a Landau level, makes use of the momentum space representation, 
Eq. (j3J), and the integral over the transverse momentum direction is thereby given as 

J dq e ± q e ± e-^ )2/2X = 2A . (9) 

From this and Eq.flS]) one finds that the factors of 7 from the Landau level electrons cancel. 

Recognizing that the nucleons are in thermal equilibrium, with Fermi momentum pp, the 
angular integrals can be done using the the pion exchange momentum in the direct diagram, 
as an independent vector by inserting: 

J d^ik-px+ps) = I . (10) 

For the D-D diagram the only angular integral needed is dQk, but for the E-E and D-E terms 
two angular integrals are needed. We use dfl^, which is almost independent, and also dfl p 
as a check. The dQ q e is modified by a factor from the standard result. 
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With no nucleoli polarization the calculation is greatly simplified, and we find 

M A _ A = 1.38 x 10* A 2 q e qzPF ( k ; kC ) formalism 
M A _ A = 3.8 x 10 4 A 2 q e q v z p% (k,p) formalism (11) 

The result for the energy integral over the Fermi distributions is that of Friman-Maxwell: 

I = 5.68 x 10 3 (A;T) 8 (12) 

Using the standard values for the nucleon Fermi momenta, and q e = 85MeV, q v = A.7kT 
we find the emissivity 

T 

e AS {k, k e ) ~ 1.18 x 10 24 (— — -) 7 erg cm^s" 1 

10 y ii 

e AS (k,p) ~ 2.96 x 10 24 (^— ) 7 erg crn^V 1 , (13) 

in the two formalisms. 

During 10s < t < 30s, 10 9 i^ < T < lO 11 -^ approximate period and temperature, the 
mean free path of the neutrinos emitted is or the order of a few cm. Thus the modified 
URCA asymmetric emission takes place when the neutrinosphere is just within the neutron 
star, as depicted. Therefore, the asymmetric energy emitted in the 10-20 s period is 






Region for modified URCA process with netrinos emitted 



E AS (k,k*) ~ 1.21 x lO^j^YiRl - Rl) ergs 
E AS (k,p) ~ 2.74 x E AS (k,k e ) , (14) 

with a factor of .52 for the neutrinos that miss the protoneutron star. 
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3.1.1 Radius of neutrinosphere during modified URCA emission 

The final step in our derivation is to estimate the volume in which the neutrino emission 
takes place with the modified URCA process in a strong magntic field, which means finding 
the radius of the neutrinosphere during the period of emission. 

Our starting point is the energy-momentum tensor, T^ u for the neutrinos with a distribu- 
tion function f v (x, k, t) for each type of neutrino, giving an energy density,?/ and momentum 
density F 

F ' = r °' = / wr k ' f - ■ (15) 

Using the Conservation of momentum, d u T 0u = gives t-dependence: 

d t U + A • F = . 

Based on these equations, and using the Spherical Eddington Model for neutrino atmosphere 
(Schindler-Shapiro, ApJ 259 (1 982) 311; Janka-Raffelt, PR D59 (1998) 023005, Barkovich- 
Olivo-Montemayor , |hep-ph /0503 1 13 ) with techniques that can be found in Barkovich et al, 
we find for the radius of the neutrinosphere during the period of emission: 

R v ~ 9.96 km (16) 

when when the temperature is in the range T ~ 10 10 A giving 

E AS ~ 1.45 x 10 42 ergs(-^r) 7 , (17) 

10 9 

and a pulsar velocity v = 1000 km/s for an expected T=10 10 K. Our results are as shown: 



lOOOOO 

v ( km/s) 

10000 
1000 



OS 1.0 1.5 

T ( 10 10 K) 

Pulsar velocity vs T assuming constant emission volume 
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4 CONCLUSIONS 



At the time the radius of the neutrinosphere is that of the neutron star the modified URCA 
process dominates the energy emission. 

At this time the electrons produced in the modified URCA process are in Landau levels, 
causing asymmetry of all neutrinos emitted. 

When the radius of the neutrinosphere is about 9.96 km vs the neutron star radius of 10 
km, the neutron star recives a kick of 1000 km and more, consistent with observation 

WE PREDICT A STRONG CORRELATION BETWEEN THE PULSAR'S VELOCITY 
AND LUMINOSCITY FOR LARGE PULSAR KICKS, AS OBSERVED (see Hansen et al) 




WE PREDICT THE BIG KICK STARTS AT ABOUT 10 s, CONSISTENT WITH AN 
ANALYSIS OF SN1987 (K. Hirata et al, P.R.D38, 448 (1988)). 




From supernova SN1987A 
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